Weinberg power counting and the quark determinant at small chemical potential 
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We construct an effective action for QCD by expanding the quark determinant in powers of the 
chemical potential at finite temperature in the case of massless quarks. To cut the infinite series 
we adopt the Weinberg power counting criteria. We compute the minimal effective action (~ p 4 ), 
expanding in the external momentum, which implies the use of the hard thermal loop approximation. 
Our main result is a gauge invariant expression for the phase 8 of the functional determinant in QCD, 
and recovers dimensional reduction in the high-temperature limit. We compute, analytically, {9 2 } in 
the range of p <C 2ttT, including perturbative and nonperturbative contributions, the latter treated 
within the mean field approximation. Implications for lattice simulations are briefly discussed. 
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I. INTRODUCTION 

There has been an increasing interest in the last few 
years in the sign problem or phase problem in QCD p|. 
For a finite chemical potential, /i, the fermion determi- 
nant matrix is non-positive definite, so it is not possible 
to perform Monte Carlo simulations in the usual fash- 
ion [3, Nevertheless, the Glasgow method [H and 
rewheighting techniques [5[ have provided great advances 
in the description of phase transitions on the lattice, con- 
sidering a set of parameters near the transition line. 

There is special interest in the region of high temper- 
ature and low chemical potential, since it corresponds to 
the sector of the phase diagram of strong interactions 
probed by high-energy heavy-ion collision experiments 
6| . In this regime it is possible to expand the fermion de- 
terminant in powers of \ij A, where A is some mass scale 
related to the temperature, A ~ T 04IH- In associa- 
tion with the ones mentioned previously, this technique 
is very convenient and successful to describe this region. 
Besides, several other complementary approaches were 
proposed with the intent to shed some light on the sign 
problem [H-GjJ]. 

In this paper we investigate the small chemical poten- 
tial sector of the QCD phase diagram. In particular, we 
present a scheme that is valid for a wide range of temper- 
atures in the soft region. Methods that are based on an 
expansion in \i always have to resort to approximations to 
compute the coefficients of the series. For this purpose, 
there are numerical approximation techniques that yield 
good results but do not allow for a deeper analytic study, 
the most freque ntly used being the random noise method 
(see, e.g., Ref. [18]). Analytic treatments are usually re- 
stricted to very high temperatures. For instance, the 
dimensional reduction effective action fl9l - [24j is valid for 
T > 2T C , where T c is the critical temperature. A com- 
plete study of the pressure within perturbative QCD for 
all values of the temperature and chemical potential can 
be found in Ref. [25| where the authors make use of 
the hard thermal loop, hard dense loop and dimensional 
reduction approximations. 



Here we start from a different perspective: instead of 
expanding all quantities and desired observables in pow- 
ers of yLt/A, the idea is to keep the relevant terms in the 
effective action according to the Weinberg power count- 
ing criteria 26] . In what follows we construct an effective 
action for QCD by expanding the quark determinant in 
powers of the chemical potential at finite temperature in 
the case of massless quarks. We compute the minimal 
effective action expanding in the external momentum up 
to order ~ p 4 in power counting. In practice, the momen- 
tum expansion performed here is equivalent to the hard 
thermal loop (HTL) approximation [27U29j . Our main 
result is a gauge invariant expression for the phase angle 
9(fi) of the functional determinant in QCD, which can be 
written as detM(/x) = |detM(/x)| e ie ^'. An interesting 
analysis of the angle 9 has been recently performed using 
the random matrix framework [ill [l7| ■ 



As a first test of our method, we recover the dimen- 
sional reduction approximation in the limit of high tem- 
perature. As a second step, we compute analytically (O 2 ) 
in the range of p <C 2ttT, keeping nonzero mode contri- 
butions in the effective action. In this calculation, we 
include perturbative and nonperturbative contributions, 
treating the latter in the mean field approximation. Im- 
plications for lattice simulations are also discussed. 



This paper is organized as follows. In Sec. II. we de- 
velop the general expansion for the fermionic determi- 
nant, present the Feynman rules and the relevant dia- 
grams, settling the framework. In Sec. III. we discuss the 
power counting hierarchy in order to be able to cut the se- 
ries and compute the minimal effective action. As a first 
test of our framework, we also reobtain the result from 
dimensional reduction in the limit of high temperature, 
including the non-zero mode terms. In Section IV we 
present the calculation of (9 2 ) in the range of p <C 2ttT. 
Section V contains our conclusions. 
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II. THE EXPANSION 

The generating functional for QCD with massless quarks 
at finite chemical potential is defined, in euclidean space, 
as 



Z 



2X?det(-0 + i/i7 4 )e~ 5YM [ G l , 



(1) 



where G are the gluon fields, also present in the covariant 
derivative = — iG^, and S YM is the Yang-Mills 
(YM) action. We can expand the fermion determinant 
in powers of the chemical potential assuming that /i < 
A ~ T: 
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where Jgdy = J Q dy A J d 3 y, S(yi>,y a ) is the dressed 
fermion propagator, which can be expressed as a series in 
powers of the gauge field and the free fermion propagator 
using the self-consistent relation 

S(x,y) = S F (x-y) - / dzS{x,z)0(z)S F {z - y) . (3) 
J0 

The expansion, then, will contribute to additional 
terms in the effective action S c!! = SVm + J2n s ^^'^ ■ 
The new terms, expressed in momentum space, are of 
the form 

= ^4dp 1 ---dp n {2T:)H{p 1 + ---+p n ) 



l ({Pi})ti^ /11 (pi).-.G Mn (p B ), (4) 



i(s,n) 



where the last integral denotes also the sum over bosonic 
Matsubara frequencies fdp = T J ,~? 3 J2 P4 =2n p irT' ano - 
5 n o- In this way, one obtains a positive- 
definite fermion determinant, and the contribution from 
the chemical potential will be part of an effective gluon 
action. 

The Fcynman rules in momentum space for calculating 
the different effective vertices are almost the same as the 
usual ones. The difference is that all of the operators be- 
tween chemical potential insertions must be transposed 
in order. Figure [T] shows a general diagram with chemical 
potential and gluon insertions. To construct a diagram 
for a vertex with n gluons and s chemical potential in- 
sertions, one puts — i ai 7Aii for any gluon insertion, —^74 
for any chemical potential insertion, and divides by the 
symmetry factor s. 

Between chemical potential insertions, the order of the 
operator must be transposed considering momentum con- 
servation. In a piece of the effective vertex shown in Fig. 
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FIG. 1: General diagram for the construction of the effective 
vertices. The small circles correspond to chemical potential 
insertions. 



FIG. 2: A set of operators between two chemical potential 
insertions. 




[21 the integrand must be written as 



S , F(fc)(-t Ql 7 AI1 )S , F(fc + pi) ■ 



{-t ar -f llr )S F {k + pi H Vp r ) 



(t) 



H74), (5) 



where the exponent (t) in the brackets is a re- 
minder to transpose the order of the operators: 

[0102 ■ ' ' Or-lO r ]^ = a0r-l • • • 0201- 

Finally, one takes the trace over gamma matrices 
and color group representation, integrating over internal 
fermionic momentum (odd Matsubara frequencies). The 
sum of all diagrams will produce the effective vertices 
r("' s ), which will be invariant under any cyclic change in 
the set of indices /ij, dj, Pi (or Xj, in the case of configu- 
ration space). 



III. THE MINIMAL EFFECTIVE ACTION 

Now, we need a criteria to cut the series. For low- 
energy effective theories, one can consider the Weinberg 
power counting p6| . which uses the argument that all 
mass parameters (external momentum, chemical poten- 
tial and gluon fields) must be less than a certain scale 
that is proportional to the temperature. This approach 
has provided a very successful description within chiral 
effective models of QCD. Although this is not exactly the 
case here, it is reasonable that for low-energy processes 
at high temperature typical values of the operators men- 
tioned above, as well as the chemical potential, can be 
considered to be in a region of the same order or smaller 
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than the scale. Then, assuming p ~ G ~ p, and expand- 
ing the effective Lagrangean in soft modes, we can cut it 
at a given order in powers of the momentum scale. 

In the case of Yang-Mills theories, this soft-mode ex- 
pansion for high temperatures corresponds to the HTL 
approximation. The minimal action must be of order 
(p) 4 . So, applying the power counting criteria, the min- 
imal effective action is given by S^ n = S VM + S^ 0,2 ' + 
5(0,4) + [ 5 (2,2) + 5.(3,1)] where ° the indices s ) are 

defined in Eq. Q, and the last two terms are expanded 
in momentum up to zeroth order which leads to the ap- 
pearance of functions of \p\/pi- 

The whole series of gauge fields is gauge invariant at 
each order in the expansion in p, i.e. ^ n S^ n ' s ^ is gauge 
invariant for all values of s as can be seen directly from 
Eq. @. Since the terms in the sum are traces contain- 
ing dressed propagators, which are gauge invariant, every 
term in the sum is gauge invariant. Moreover, the mini- 
mal effective vertices that we need satisfy Ward identities 
of the form 



(6) 



(7) 



and analogous relations obtained by changing cyclical 
indices and arguments. Equation (J7]) vanishes, since 
pC 2 * 1 ) — 0. As is well known, HTL preserves the Ward 
identities. 

The nonvanishing diagrams for the vacuum contribu- 
tions r( s, °) are the known vacuum corrections to the ther- 
modynamic potential 



r(°' 2 ) 
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For s > 4 all contributions vanish, as was demonstrated 
in Ref. (30j . The next nonvanishing term has the form 
of the polarization tensor in the HTL approximation 
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with the lightlike four- vector k — (k,i). Finally, the 
vertex components which correspond to 16 are 
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(10) 



The dimensional reduction approximation can be ob- 
tained directly from the last expressions of the effec- 
tive vertices by simply considering the case in which 
|p| <C 2ttT. In order to expand in powers of the external 
momentum we have to separate the zero mode from the 
other modes in the gluon fields G p = Ap + Bp, with 



A p (x) = T 



d 3 p 



{2n)3 e^Gp(p,0), (11) 
Bp(x) = T ]T J ' ^ e*" Gp{p), (12) 



Ap being the field corresponding to the zero mode. The 
lowest order contribution in the expansion {p — > 0) is 



given by 



S< 2 



5(3,1) = 



2ir 2 
3tt 2 



dx tr 



Ai - ~B 2 







dx tr [A\ + A 4 B 2 ] 



(13) 



(14) 



If we set B = 0, we recover the dimensional reduc- 
tion effective action at tree level 0, [24j]. The usual 
effective action for the high-temperature regime is con- 
structed through loops integrating the B fields. However, 
this is done perturbatively, and this is not the regime in 
which we are interested. The fact of including B as a 
nonperturbative field (except for the power counting) in 
principle will enhance the range of validity in tempera- 
ture to values lower than in the case of just considering an 
infinite temperature expansion. In this sense, our frame- 
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work goes beyond dimensional reduction, and can probe 
temperatures closer to T c . 



IV. 



2 ) AT HIGH TEMPERATURE 



The phase angle 9 of the complex functional determi- 
nant is a crucial quantity for lattice simulations in QCD. 
The knowledge of 9 allows for the separation of the func- 
tional integral into two different regions: 16*1 tt/2. Nev- 
ertheless this is not a simple task, since it is not restricted 
to values of — tt < 9 < it but seems to increase with the 
volume [3l| . The calculation of the average phase 



2i8\ 



det(— ilp + i^-fi) 
det(- 



(15) 



gives a measure of how problematic could be the phase in 
lattice simulations. As this quantity must vanish in the 
thermodynamic limit, it should be evaluated in a finite 
volume (otherwise the average phase above vanishes as 
soon as /x is nonzero). The random matrix aproximation 
QJ] yields a vanishing result for fj, > m^/2. One can 
use our effective expression for the phase to compute the 
average phase factor in Eq. (fl~5"j) . If one expands the ex- 
ponential in powers of the angle, the relevant term will 
be ~ {9 2 ), since the average of the phase must be real. 
In particular, the determination of (9 2 }p, the average of 
the angle as a function of the plaquette, is an impor- 
tant ingredient in order to localize the critical line in the 
temperature-number density phase diagram of QCD [3lj |. 

In our framework, we have S^ 3 ' 1 ) = id. Writing the 
fields in the form 



A«(x) = T f dx^x) 
Jo 



(16) 



Bftx) = G°(x)-T / dx A Gl{x) , (17) 
Jo 

we can use Eq. (1141) to express 9 in the high-temperature 
limit as 

P 



J^I' 1 "'" ' I << ''■'■ J 'I'' .'hi ;(■■'■ (y)G b 3 (x if,".; ( r i 



(18) 



where x = (a;, £4), y = (x, 2/4) and d abc are real and 
totally symmetric as usual. The expectation value of 9 
vanishes, so we compute the next power which can be 
expressed in terms of two-point correlation functions as 



fi 2 T 2 N 2 (N 2 -4) r? 
(12tt 2 ) 2 N c (N 2 - l) 2 J 



dx^dx\dy^dy\ 



x / d 3 xd 3 x'{2(Gt(y)Gt(y')}(Gt(x)Gt(x')) 2 
Jv 

+4(GUy)GUx'))(Gt(x)GUy'))(Gt(x)GUx'))} (19) 

where the expression is not exact since we allow for non- 
perturbative contributions, and we used (G a G h ) ~ S ab . 



If we separate the gluon field in a perturbative and 
a nonperturbative contribution, G = G p + G np , with 
(G p G np ) — 0, there will be contributions to (9 2 ) com- 
ing from a purely perturbative term, a purely nonpertur- 
bative term and the crossed terms. The basic building 
blocks are the two-point functions for the field strengths, 
which we compute in the sequel. 

For the perturbative case, we start with the HTL ef- 
fective Lagrangian 



+—C a 



1 



TjidtiG®) 2 



Ga 
R, 



Ait Ui . n\2 ' 



(ik ■ D) 

where m is the standard HTL effective mass [3 
^ 2 r ZN f ^ 2 
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T\N f + 2N C ) + 



(20) 



(21) 



Then, the fourth component of the gluon two-point func- 
tion can be written, in the limit p <§; 2ttT, as 



{G%{x) p G\{x') p ) 



5 ab g 2 T 
1 



ip.(x— as') 



E 

n#0 



d 3 p 
J2nf" 



(22) 



where uj n — 2nirT is the bosonic Matsubara frequency. 

For the nonperturbative case it is convenient to use the 
Schwinger gauge, also known as fixed point or coordinate 
gauge [331 ] . x ■ G = 0. Then, the gauge field can be 
expressed in terms of the field strength tensor: 



G a Jx) 



dssG® a (sx)x c 



(23) 



If we approximate the two-point function for the field 
strength by its mean field value, i.e. 



(GUx)GUy)) » (GL(O)GL(O)) = ~{£ 2 ) 
where £f is the color electric field, we find 

{G2(x) np G1(x') np ) ~ ~Y2~ X ' X ' ' 



(24) 



(25) 



Using the results above, the computation of (9 2 ) is long 
but straightforward. Collecting all terms, and assuming 
a large volume, we find 



^ 2 N 2 (N 2 ~4) 
(12tt 2 ) 2 N c 
3(£ 2 )ff 4 R 5 
32ir 5m 



TT(£ 2 ) 2 g 2 fW_ R? 
8(N 2 - 1)T \m 2 + m 4 

+ F(/3™,/3A M g)T 5 6 i? 3 



(26) 



where R is the radius of the system in spherical coordi- 
nates, A^jg is the energy scale in the modified minimal 
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subtraction (MS) scheme, and F is an integral that can 
be computed numerically. 

The last term is purely perturbative, whereas the other 
ones come from the mixed contribution. The purely 
nonperturbative contribution vanishes identically in the 
mean field approximation. Notice that, in the thermo- 
dynamic limit, the dominant contribution is the one pro- 
portional to R 7 ~ V 7 / 3 , so that (9 2 ) grows faster than 
quadratically with the volume, the other contributions 
being essentially irrelevant. 

V. CONCLUSIONS AND OUTLOOK 

We presented a well-defined and systematic procedure 
for computing the fermionic determinant in QCD at fi- 
nite temperature and chemical potential. In the frame- 
work defined by the power counting method described 
above, we calculated exactly the minimal corrections to 
the effective action. This procedure can be systemati- 
cally extended to higher orders. 

We also computed explicitly (9 2 ) in the high- 
temperature limit and analyzed its volume dependence. 
In this limit, the dominant contribution comes naturally 
from the zero Matsubara mode, so that dimensional re- 
duction is a good approximation. Although we do not 
present results beyond the mean field approximation for 



the gluon condensate, we expect that in that case, a new 
length scale associated with the nonuniformity of the con- 
densate will compete with the scale R but should not 
modify the picture appreciably in the limit considered in 
this work. 

In this paper we proposed a new way of approach- 
ing the problem of determining the phase of the complex 
functional determinant in QCD at small chemical po- 
tential and relatively large temperatures, an important 
issue for lattice QCD thermodynamics. Our framework 
clearly contains, and goes beyond, the dimensional reduc- 
tion approximation, being valid for temperatures closer 
to T c . We believe that this analysis, complemented by 
previous existing results, might shed some light onto the 
phase transition in the region of high temperature and 
low density presumably probed by relativistic heavy-ion 
experiments, and perhaps indicate a simple way to han- 
dle the sign problem for temperature values closer to the 
critical temperature. 
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